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The magnetic analog of the Griineisen parameter, i.e., 
the magnetocaloric effect, is a valuable tool for studying 
field-tuned quantum phase transitions. We determine the 
magnetic Gruneisen parameter of the one-dimensional 
random transverse-field Ising model, focusing on its low- 
temperature behavior at the exotic infinite-randomness 
quantum critical point and in the associated quantum 
Griffiths phases. 



We present extensive numerical simulations showing 
that the magnetic Griineisen parameter diverges loga- 
rithmically with decreasing temperature in the quantum 
Griffiths phase. It changes sign right at criticality. These 
results confirm a recent strong-disorder renormalization 
group theory. We also compare our findings to the be- 
havior of the clean transverse-field Ising chain. 
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1 Introduction Zero-temperature quantum phase 
transitions are characterized by a competition between 
different quantum ground states. The resulting quantum 
fluctuations often cause unconventional behavior of ther- 
modynamic and transport properties in the vicinity of the 
phase transition. Celebrated consequences include non- 
Fermi liquid behavior in metals and the emergence of new 
phases such as exotic superconductors (for recent reviews 
see, e.g., Refs. dEHUBl). 

Realistic systems often feature considerable amounts 
of quenched randomness in the form of vacancies, impuri- 
ties, or other defects. In these systems, the combination of 
static random fluctuations and quantum fluctuations leads 
to dramatic disorder effects in the vicinity of a quantum 
phase transition. They include quantum Griffiths singulari- 
ties |l5]|6|, infinite-randomness critical points with activated 
dynamical scaling |I7][8l, and smeared quantum phase tran- 
sitions ll9l fT0l . A recent review of these phenomena can be 
found in Ref. ifTTl . 

In recent years, the Gruneisen parameter F 01211131 
has been established as a valuable tool for characterizing 
pressure-tuned quantum phase transitions. It is defined as 
the ratio between the thermal volume expansion coefficient 
P and the specific heat Cp. If the quantum phase transition 



is tuned by varying an external magnetic field, one needs to 
consider the magnetic analog Fh of the Gruneisen parame- 
ter, defined as the ratio between the temperature-derivative 
of the magnetization and the specific heat ch- It can be 
directly measured via the magnetocaloric effect. 

According to Zhu et al. 01311141 . the thermal expan- 
sion coefficient is more singular than the specific heat at 
a generic pressure-tuned quantum critical point, leading to 
a diverenge of the Gruneisen parameter. In particular, these 
authors showed that F ^ if the temperature T is 

lowered at the critical pressure pc and F ~ l/(p — pc) 
if the pressure p approaches its critical value at zero tem- 
perature. Here, z and i' are the dynamical and correlation 
length critical exponents, respectively. These results hold 
below the upper critical dimension (where hyperscaling is 
valid). Above the upper critical dimension, F still diverges, 
but the functional form depends on dangerously irrelevant 
variables. The divergence of the Gruneisen parameter has 
been confirmed experimentally at several pressure-tuned 
quantum phase transitions, and the functional form of the 
divergence can be used to discriminate between different 
theories QsKHKIIKIII . 

The above results apply to clean quantum critical 
points with conventional power-law scaling. Given the 
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utility of the Gruneisen parameter, it is very desirable 
to know its behavior within the unconventional quantum 
phase transition scenarios emerging in the presence of 
quenched randomness. Recently, the low-temperature be- 
havior of Griieneisen parameter F and its magnetic analog 
Fh were studied at infinite-randomness quantum critical 
points and in the associated quantum Griffiths phases |fT9l . 
Using heuristic rare region arguments as well as a scal- 
ing theory, it was predicted that the Gruneisen parameter 
diverges logarithmically with decreasing temperature T 
in the quantum Griffiths phases and as a higher power of 
ln(l/T) right at the infinite-randomness quantum critical 
point. 

In this paper, we confirm and illustrate these predic- 
tions by considering a specific microscopic model featur- 
ing an infinite-randomness quantum critical point, viz., the 
random transverse-field Ising chain. Since the transition 
is tuned by magnetic field, we study the magnetic analog 
Fh of the Gruneisen parameter (i.e., the magnetocaloric 
effect). Our paper is organized as follows. In Sec. |2l we 
define the model, and in Sec. [3j we briefly summarize the 
results of the scaling theory of Ref. |fT9l . Section |4] is de- 
voted to extensive numerical calculations of the magnetic 
Gruneisen parameter. We conclude in Sec. |5] 

2 Random transverse-field Ising chain We con- 
sider a one-dimensional random transverse-field Ising 
model given by the Hamiltonian 



H = 



La 



(1) 



Here, erf and erf are the Pauli matrices representing the 
X and z components of the spin- 1/2 at site i, and Lq is 
the chain length. The bonds Ji > and transverse fields 
hi > are independent random variables with bare (initial) 
probability distributions Pi (J) and Rj{h). The physics 
of this model has been studied extensively in the litera- 
ture. If the bonds are sufficiently strong compared to the 
fields, the ground state is ferromagnetic, with a sponta- 
neous magnetization in ±z direction. In the opposite limit, 
for sufficiently strong fields, the ground state is param- 
agnetic with vanishing ^-magnetization. The two ground- 
state phases are separated by a quantum critical point lo- 
cated at [ln( Ji)]dis = [ln(/ii)]dis where [. . .]dis denotes the 
disorder average. The model does not show long-range or- 
der at any nonzero temperature. 

The quantum critical behavior can be determined ex- 
actly by means of the strong-disorder renormalization 
group llllSl. It does not depend on the details of the dis- 
tributions Pi (J) and Ri{h) provided they are not too 
singular The quantum critical fixed point that emerges 
from this analysis is of exotic infinite-randomness type. 
This implies that correlation time and correlation length 
^ are related by the exponential law In(^T-) ^ S,"^ with 
a tunneling exponent i/; = 1/2. In contrast, the relation 



between the correlation length and the distance r from 
criticality is conventional, ^ ~ kl with a correlation 
length exponent v = 2. 

In the following, we consider systems with randomness 
in the bonds Ji only. The transverse magnetic field is taken 
to be uniform, Ri{h) = 5{h — H) because we wish to view 
it as the experimental tuning parameter. 

3 Magnetic Gruneisen parameter The magnetic 
analog Fh of the Gruneisen parameter is defined as the 
ratio of the temperature derivative of the magnetization 
M = J^ii'^i) ^^'^ specific heat, both taken at con- 
stant field H IfTSl . Note that M is the x-magnetization (the 
magnetization conjugate to the external field) rather than 
the ^-magnetization (the order parameter of the transition). 
Because the order parameter vanishes at any nonzero tem- 
perature, M is the magnetization measured in experiment. 
We thus define 



F, 



H 



idM/dF)H 



c 



H 



{dS/dH)T 
T{dS/dF)H 



1 

T 



dT 
dH 



s 
(2) 

with S the entropy. Here, {dS/dH)T explores the depen- 
dence of entropy on the dimensionless distance from criti- 
cality, r = [H — i/c)/^^c-ElThe last equality in ^ shows 
that Fh can be obtained directly from the magnetocaloric 
effect. 

Let us now apply the scaling theory |fT9l for the mag- 
netic Gruneisen parameter to our case. According to the 
strong-disorder renormalization group fT^, the entropy is 
proportional to the number of independent active clusters 
at energy scale F. It thus has the scaling form (see, e.g., 
Ref. lEUl) 

S{r,F) = LoSo[ln(To/T)]-2<^(r ln(To/T)) (3) 

with So = ln(2) the entropy (in units of ks) of a single 
cluster and Fq a microscopic temperature scale. The scal- 
ing function 'P{y) is analytic because there is no phase tran- 
sition at any finite temperatures. For small y, we therefore 
expand <P{y) = <P{0) + y<P'{0) + y^<P"{0) + .... At the 
quantum critical point, the random transverse-field Ising 
model is self-dual (invariant under the exchange hi ^ Ji). 
The scaling function is therefore even in y which im- 
plies <P'{0) = 0. For large y, i.e., in the Griffiths phase, 
^(y) = Ay^ exp{—c\y\) with A and c constants. 

We first consider the quantum critical region defined 
by |r I ln(To/T) ^ 1. Taking the appropriate derivatives of 
the entropy and keeping only the leading terms, we find the 
specific heat 



CH = 2LoSo<PmHFo/F)]- 



(4) 



' Within the strong-disorder renormalization group (7][8l, the 
distance from criticality is naturally measured in terms of 5 ~ 
[ln(/ii)]dis — [lii(>/i)ldis- Close to criticality, r and S are propor- 
tional to each other. Using S instead of r would thus introduce 
some constant factors at intermediate stages of the calculations, 
but leave the final results unchanged. 
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and the magnetization derivative 

{dM/dT)H^LoSo r$"{0)/H,. 
This results in a magnetic Griineisen parameter 



2m 



MTo/T)] 



(5) 



(6) 



which has a sign change right at criticality r = 0. Note 
that Fh takes a different form than in the generic theory of 
Ref. ||T9| because the first derivative of the scaHng function 
vanishes in our system. 

In the quantum Griffiths phase, \r \ ln(To/r) ^ 1, we 
use the large-argument limit of the scaling function 'P{y) 
to obtain the entropy S = Log{r){T/To)^^''^ with g{r) = 
Asor'^. The Griffiths exponent is given by A(r) = c\r\. 
Taking the appropriate derivatives leads to the specific heat 



Cff =io.9(r)c|r|(T/To)^M 



(7) 



and the magnetization derivative 



{dM/dT)H = -^9{r) ^(r/ro)^('-) ln(ro/T) . (8) 



Forming the ratio results in the magnetic Griineisen param- 
eter 

rn^j^^HnlT). (9) 

The Griineisen parameter thus diverges logarithmically 
with decreasing temperature in the entire quantum Grif- 
fiths phase. The prefactor diverges itself as 1/{H — He) 
when approaching criticality. 

4 Computer simulations In this section, we report 
the results of extensive computer simulations of the ran- 
dom transverse-field Ising chain ([T]i by means of a nu- 
merical implementation of the strong-disorder renormal- 
ization group. We have studied long chains of up to 10^ 
sites with periodic boundary conditions. All results are av- 
erages over a large number {N = 5000 . . . 25000) of dis- 
order realizations. The (bare) system is moderately disor- 
dered; the bonds Ji are drawn from the interval [0.5, 1] 
with a flat distribution. The transverse field is uniform, 
hi = H. The critical field can be determined analytically 
from the condition [ln(Ji)]dis = [ln(/ii)]dis- This leads to 
He = exp[ln(2) - 1] « 0.73576. 

We first consider the behavior of the entropy as a 
function of temperature T and distance to criticality r. 
To find the isentropes (curves of constant entropy), we 
run the renormalization group until a certain fixed num- 
ber S/ ln(2) of clusters is left. The corresponding temper- 
ature is then given by the energy cutoff of the renormaliza- 
tion group. Figure [T] shows the resulting isentropes in the 
temperature-field plane. The entropy strongly accumulates 
near the critical field. The derivative {dS/dH)T is thus 
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Figure 1 Isentropes of the random transverse-field Ising 
model ([T]). Note the symmetry of the cures w.r.t. r ^ — r. 
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Figure 2 Scaling plot of the isentropes of Fig. [T] according 
to the scaling form ^ of the entropy. 



positive on the magnetic side of the transition (H < He) 
and negative on the nonmagnetic side (H > He). This is 
responsible for the above-mentioned sign change of Fh at 
criticality. To verify the scaling form (|3]l of the entropy, we 
redraw the same data in a scaling plot in Fig. [2] The figure 
shows that the high-entropy curves (which correspond to 
high temperatures and the early stages of the renormaliza- 
tion group) do not scale perfectly and also reflect the H-J 
asymmetry of the bare system. However, with decreasing 
entropy (i.e., decreasing temperature), the scaling becomes 
better and better until it is almost perfect within our error 
bars. 

We now turn to the magnetic Griineisen parameter Fh- 
Figure [3] gives an overview of Fh as a function of trans- 
verse field for various temperatures. In agreement with the 
theory of Sec. |3] the Griineisen parameter is negative for 
H < He and positive for H > He- At the critical field 
He, it vanishes for all temperatures. The figures also shows 
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Figure 3 Overview of the magnetic Gruneisen parameter 
of the random transverse-field Ising model. 
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Figure 4 Magnetic Gruneisen parameter as a function of 
ln(ro/T) for various values of the transverse field. 
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Figure 5 Prefactors 7 = rH/ln(To/r) of the low- 
temperature Gruneisen parameter as a function of trans- 
verse field H. The solid line represents the theoretical pre- 
diction from (|9]l without adjustable parameters. 
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that, at a given field H, the magnitude of the Gruneisen pa- 
rameter increases with decreasing temperature. 

To investigate more closely the temperature depen- 
dence of the Gruneisen parameter in the quantum Grif- 
fiths phase, we plot Fh as a function of ln(r/To) in 
Fig. |4] In the systems furthest away from criticality 
(r = 0.00508, -0.00511), Fh follows (after a very short 
high-temperature transient) the ln(ro/r) behavior pre- 
dicted in ^ for the Griffiths phase. If we move closer to 
criticality (r = 0.00304, -0.00307), the high-temperature 
behavior of Fh is clearly faster than ln(r/ro)- At 
lower temperatures, it crosses over to the Griffiths phase 
form ln(To/T). Finally, the systems closest to criticality 
(r = 0.00101, -0.00103) do not reach the crossover to the 
Griffiths phase form within the temperature range covered 
in the simulations. 

To confirm (|9|l quantitatively, we have performed anal- 
ogous calculations for many more values of the transverse 
field and extracted the prefactors of the ln(ro/r) low- 



Figure 6 Derivative xth = [dFn / dr)T of the Gruneisen 
parameter taken at criticality, r = 0. The soHd line is a fit 
to the expected [ln(To/r)]'^ behavior at low temperatures. 



temperature behavior of the Gruneisen parameter These 
prefactors 7 = Fh/ ln{To/F) are plotted in Fig. |5] The 
overview given in the inset shows that 7 increases and ap- 
pears to diverge close to criticality. For a detailed com- 
parison with (|9]l, we plot 7^^ as a function of H — He 
in the main panel of the figure. The numerical data are 
in very good quantitative agreement with the prediction 

7-1 = H - He. 

Finally, we turn to the behavior of the magnetic 
Gruneisen parameter in the quantum critical region, 
|r|ln(ro/r) <C 1. To verify (|6]l, we compute the 

Gruneisen "susceptibility" x/^H = {dFH/dr)T at critical- 
ity. According to this quantity is predicted to behave as 
[ln(ro/2^)]"^ at low temperatures. The data in Fig. |6] show 
that xth indeed follows the expected behavior over many 
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orders of magnitude in T. The small deviations at the low- 
est temperatures can be attributed to finite-size effects as 
well as systematic errors introduced by the finite field step 
5H \n our simulations. 

5 Conclusions In summary, we have performed ex- 
tensive computer simulations of the one-dimensional ran- 
dom transverse-field Ising model by means of a numeri- 
cal strong-disorder renormalization group. We have stud- 
ied the magnetic analog of the Griineisen parameter (which 
is related to the magnetocaloric effect), focusing on the 
low-temperature behavior at the infinite-randomness quan- 
tum critical point and in the associated quantum Griffiths 
phases. 

Our results confirm the predictions of a recent scal- 
ing theory (T9l of the Griineisen parameter and its mag- 
netic analog. Specifically, we have found that the magnetic 
Griineisen parameter diverges as ln(ro /T) with decreasing 
temperature in the quantum Griffiths phases on both sides 
of the transition. The amplitude of the logarithmic temper- 
ature dependence is universally given by 1/{H — He) and 
does not depend on details of the Hamiltonian. The behav- 
ior right at critic ality differs from the generic scaling re- 
sult of Ref. lfT9l . This is caused by the self-duality of the 
random transverse-field Ising model at its quantum criti- 
cal point which results in a scaling function that is even 
in the distance r from criticality. Therefore, the leading 
scaling contribution to the Griineisen parameter vanishes 
right at criticality. The next-to-leading term produces a 
sign change with a slope that diverges as [ln(To/T)]'^ with 
decreasing temperature. In an experimental field sweep 
across the quantum critical point, this rapid sign-change 
leads a strong feature in the magnetocaloric effect. 

Let us compare the behavior of our random transverse- 
field Ising chain to that of a clean chain |T4|. The magnetic 
Griineisen parameter of the clean chain does not diverge 
with decreasing temperature but saturates at sufficiently 
low T (even though the saturation value can become ar- 
bitrarily large close to criticality). The critical point of the 
clean chain has the same self-duality as that of the random 
chain. Therefore, analogously to our case, the leading scal- 
ing contribution to the Griineisen parameter (which would 
be proportional to T^^) vanishes in the quantum critical 
region. This results in a sign change right at the critical 
field with a Griineisen "susceptibility" xfh ~ T~^. 

The random transverse-field Ising chain studied here 
provides one of the simplest examples of a field-tuned 
quantum phase transition in the presence of disorder It 
is important to note that the model realizes random- Tc 
(random mass) disorder This means the disorder does not 
(locally) break the order parameter symmetry, it just fa- 
vors one phase over the other. In many realistic disor- 
dered systems, the interplay of an external (uniform) mag- 
netic field and impurities and defects can generate stronger 
random-field type effects 121112211 . In particular, this ap- 
plies to LiHo2:Yi_2;F4 which is reasonably well approxi- 



mated by a transverse-field Ising model otherwise 123112411 . 
If strong enough, such random fields invalidate the infinite- 
randomness critical point scenario and require a separate 
analysis. 
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